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Abstract 

Suppose  that  Xz-j^x  <  —  <  X/  •>  are  the  order  statistics 
of  a  sample  of  size  n  from  a  continuous  and  symmetric  distri- 
bution with  a  location  parameter  9. 

We  consider  a  (fictitious)  subsample  of  size  k  drawn 
randomly  from  the  order  statistics  of  size  n^  and  let 
Y/-,\  <  ...  <  Y/,  \  be  the  order  statistics  obtained  from  the 
subsample. 

We  shall  determine  the  coefficients  of  a  linear 
combination 

Y~   c   Y/  -. 

^ a  [a) 

so  as  to  minimize  V(  5~  c  Y ,    s    \    X, ^  ^    <    .  .  .    <  X,    \)   imder 

^  '^ —  a   (a)  '   (1 J  (n) ' 

the  condition  that  y~  c      =   1.      Then  we  take  the  conditional 

' —  a 

expectation  and  regard  it  as  an  estimator  of  9, 
r  =  E(  ^c^  Y(^^  1  X(,)  <  ...  <X(^)) 


=  >   c   T 

^ a  a 

where,  actually,  T^  =  E(Y/^J  X/-^x  <  ...  <  X/  x)  is  a  linear 
combination  of  the  X/^%'s.   For  a  wide  class  of  distributions 
satisfying  a  certain  set  of  regularity  conditions  it  is  shown 

A 

that  9   is  asymptotically  efficient,  when  n  — >  oo  and 
simultaneously,  k  — >  oo  in  an  appropriate  way. 

Monte  Carlo  experimentation  shows  that  this  estimator, 
or  rather  a  modified  version  with  symmetric  coefficients  i.e. 
c^  =  ^]^_^^i>   attains  high  relative  efficiency   for  several 
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varieties  of  distrib  .       ven  for  small  samples  of  size 
n  =  10,  15,  20  or 

We  can  also  find  a  natural  way  of  "Studentization" 
and  asymptotic  confidence  intervals  are  obtained.   Utilizing 
Monte  Carlo  techniques,  these  intervals  also  seem  to  be 
satisfactory  for  small  samples. 
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§1.   Introduction. 

Most  of  the  so-called  "robust"  estimators  for  location 
parameters  are  either  linear  Gombinations  of  order  statistics 
or  closely  related  to  them.   Trimmed  and  Winsorized  means 
(Tukey  [28])  are  the  most  extensively  studied  examples 
(Bickel  [h],    Gastwirth  and  Cohen  [12]^,  Grow  and  Siddiqui  [10]). 
Also,  minimax  linear  unbiased  estimators  have  been  discussed 
(Birnbaum  et  al.  [5,6],  Gastwirth  and  Rubin  [13 J^  Yhap  [29]). 
Ruber's  estimator  (17]  is  primarily  of  another  kind,  but  is 
reduced  to  a  kind  of  trimmed  mean   for  his  specific  choice 
of  his  p  function.   Hodges-Lehmann' s  [l6]  is  also  of  another 
kind,  and  seems  to  be  the  only  practically  tractable  estimator 
derived  from  the  rank  tests.   (For  general  estimators  based 
on  rank  tests  see  Adichie  [1].)   But  it  has  also  been  shown 
that  it  is  also  closely  related  to  a   linear  combination  of 
order  statistics. 

It  has  been  known  since  Bennett  [2],  Jung  [19]  and 
Blom  [7]  that  the  best  linear  unbiased  estimator  (BLUE)  is 
asymptotically  efficient  in  the  sense  that  the  variance  of 
the  BLUE  is  asymptotically  equal  to  the  Cramer-Rao  bound 
for  the  variances  of  the  unbiased  estimators  (Jung  [19])j 
or  asymptotic  variance  of  the  BLUE  is  equal  to  it  (e.g. 
Chei'noff  et  al.  [8]).   Of  course,  in  order  to  calculate  the 
coefficient  of  BLUE,  it  is  necessary  to  compute  the  variance- 
covariance  matrix  of  order  statistics,  but  Blom  [7  ]  also 
showed  that  asymptotically  the  coefficients  c(iln)  of  the 
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i-th  ■ Istic   in  a   sar:  if  size  n. 


,  ,  (1.1) 


c(i  n)  oC  -        f(x) 


•«..crc  i  a-.-..-.'-.  L..-:  a-:i.^^y  and  i{    )  the  quantilc  ^i    •.,.- 
distribution.   This     ; ossion  has  suggested  that  if  we 
have  a  very  larg        ,  we  might  be  able  to  estimate 
the  density  frunction  from  the  sample,  and  then  construct 
an  estimator  which  is  asymptotically  equivalent  to  the  BLUE 
without  assuming  knowledge  of  the  shape  of  the  distribution. 
Stein  [22]  and  Hajek  [I5]  investigated  such  possibilities 
(also  Bhattacharya  [3  ]).   However,  since  in  the  above 
formula  the  second  order  derivative  of  the  density  is 
involved,  it  is  intuitively  obvious  that  the  sample  size 
must  be  prohibitively  large  in  order  that  the  asymptotic 
efficiency  be  attained. 

There  are  many  other  possible  ways  of  using  the 
information  contained  in  the  sample  configuration  to  choose 
a  linear  estimator  so  as  to  attain  high  or  full  asymptotic 
efficiency  for  a  wide  class  of  shapes  of  the  distribution. 
Such  estimators  may  be  called  quasilinear  estimators; 
one  example  was  given  by  Hogg  [I7]. 

In  this  paper  we  shall  introduce  a  quasilinear 
estimator  which  is  shown  to  be  asymptotically  (nearly) 
efficient  for  a  wide  class  of  symmetric  distributions  with 
certain  regularity  conditions.   Monte  Carlo  experiments 
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show  that  it  has  quite  high  relative  efficiencies  even  for 
moderate  sample  sizes  such  as  n  =  10,  1^,    20.   Further, 
there  is  a  natural  method  of  "Studentization"  for  this 
estimator,  which  is  shown  to  be  asymptotically  valid  for  any 
(regular)  shapes  of  distributions,  and  Monte  Carlo  results 
also  indicate  that  for  small  samples  the  resulting  confidence 
intervals  are  uniformly  approximately  valid,  and  more  or  less 
biased  on  the  conservative  side  when  the  underlying  distri- 
bution is  long-tailed. 

The  basic  idea  of  the  method  is  as  follows. 
Although  it  is  quite  impossible  to  estimate  the  variance- 
covariance  matrix  of  the  order  statistics  directly  from 
the  sample,  it  is  possible  to  estimate  the  variance-covariance 
of  the  order  statistics  of  a  sample  of  size  k,  from  a  sample 
of  size  n,  and  if  n  is  appreciably  larger  than  k,  the 
estimates  would  be  accurate.   But  if  the  distribution  is 
regular,  this  information  about  the  variance-covariance 
matrix  of  the  order  statistics  of  a  sample  of  size  k  can  be 
effectively  used  to  construct  a  "nearly"  efficient  estimator 
based  on  the  sample  of  size  n  provided  that  k  is  also  large. 

This  method  can  be  regarded  as  a  refinement  of 
the  method  proposed  by  the  author  [24]  which  was  based  on 
random  division  of  the  samples  into  several  sets  of 
subsamples,  and  also  asymptotically  uniformly  efficient, 
but  which  had  an  undesirable  property  of  not  being  a  function 
of  the  order  statistics. 
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that  Xj^  —  X^  are  independently  and  identically 
dlst  .dom  varl         .   atinuous  distribution  with 

a  ■  ■  "nmeter  0,      We  shall  denote  the  density  function 

by  i\x-d),  and  we  shall  assume  the  following: 

(1)  f(x)  >  0  for  all  X 

(2)  f(-x)  =  f(x),   that  is  the  distribution  is  symmetric 

(3)  f(x)   is  continuous  and  almost  everywhere  differentiable, 
and 

If  =  I^^T^     d-  <     ^  (2.1) 

Then,  imder  some  well  known  regularity  condi- 
tions, it  is  known  that  (Cf.  Cramer  [9])  for  any  unbiased 

A 

estimator  9   of  9, 

V(9)   >  -^        .  (2.2) 

n  I^ 

Now  let  X,-,s    <  ...  <  X/  N   be  the  order  statistics 

obtained  by  rearranging  X,  . . .  X  .   Consider  a  linear 

combination 

A    n 


e  =|Z:c,  X(,)  .  (2.3) 


-i  =  ^n-i+1   ^°^  ^  "  1.2,...   and  H  c^  =  1 


i=l 

If  c^  =  c^ 

in  (2.3),  9   is  unbiased  for  9   for  any  symmetric  distribution, 
and  we  shall  call  it  a  linear  linbiased  estimator. 

For  any  fixed  symmetric  f,  let  c^  be  defined  as 
the  value  of  the  coefficient  minimizing 

■(i) 
-h- 


V(0)  =  V(  ^  c.  x,.0  , 


under  the  condition  that  2_  c.  =  1.   Then  it  is  well  known 
(Lloyd  [20])  that  c.  is  given  by 

i   J 
where  {a^^}    =    io^.]'^   and  a^.    =   E(X^.jX^jj)  -  E(X^^^  )E(X^  ^  ^ ) 

■X-  * 

It  is  easily  shown  that  c   .  ,-,  =  c.  ,  hence 

0^  =  X:  c.  X^.^  ^         (2.5) 

is  a  linear  unbiased  estimator^  and  is  of  minimum  variance 
when  f  is  the  time  distribution,  and  is  called  the  best 
linear  unbiased  estimator  (for  short  BLUE). 

It  was  proved  by  Bennett  [2],    Jung  [19]  and  Blom  [7] 
that  under  some  regularity  conditions  BLUE  is  asymptotically 
efficient,  in  the  sense  that 

nV(g*)  — >  ^  (2.6) 

More  precisely,  for  some  class  ^   of  symmetric 
distributions  satisfying  some  regularity  conditions,  the 
convergence  in  (2.6)  is  uniform,  that  is 

nV(^*)  =  ^  +  O(^)  (2.7) 

or  still  more  precisely,  for  some  K  >  0 

nV(^*)  ^  ^  +  f     for  all  f  in  5^  .        (2.8) 

Hence  for  any  e  >  0,  if  we  take  k  large  enough,  we  have 
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^^(K)   1  T"  +  £     for  all  f  in  ^  .       (2.9) 
f 

Now  suppose  that  k  is  fi,  _      Jn.   to  satisfy 

(2.9).   •■   •  ^^  '   ^te  the  BLUE  for  th     ^    "   '    k 
specifically  by 

Now  consider  the  order  statistic  of  size  n,  n  >  k  , 
X/-,\  <  ...  <  X/  \,    and  suppose  (hypothetically )  that  a 
subsample  of  size  k  be  randomly  chosen  from  this  sample j 
the  order  statistic  obtained  from  which  is  denoted  by 

Y(l)  ^  •••  ^  "^{K)'      '^^^ 

■^(a)  ^  -^(1  )  '   a=  1,...  k,  i-j^  <  i2<  ..  .  <  ij^   (2.11) 

The  joint  distribution  of  Y,^\    <  ...  Y/,  n  is  exactly 
the  same  as  the  joint  distribution  of  the  order  statistics 
of  size  k  from  the  original  distribution.   Hence  if  we  put 

e„  =  S~~   c  Y/  ^  (2.12) 

f   fri     a   (a)  ^    ' 

we  have  kV(0^)  =  kV(^*  )  1  —  +  e  • 
I        I    i^ 

Now  we  shall  denote  the  conditional  expectation  of  6 


given  the  order  statistics  X/-,\  <  ...  <  X/  \ 


by 


.^  .      J^  * 


«?  =  E(^,  I  °n)  =  ^  ^l  \  (2-^5) 


where  T  =  E(Y/  \  0_),  and  0_  is  the  abbreviation  denoting 
a      (a) '  n'       n 

the  order  statistic  of  size  n.   T  can  be  expressed  as 
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T^  =  J^  Pa  X(,)  (2.14) 

and 

^a  =   ^r^^(a)  =  ^(i)^  "  i-l^a-1  n-A-a   ^  n\ 

for  a<_i<_n-k+a     (2.I5) 
-     0     otherwise 

Here  we  need  the  following  lemma. 

LEMMA.   Suppose  that  *(x-,  ...  x,  )  is  a  real  valued  function 

of  k  real  variables.   Let  0  =  ^-^Ci  ")  '^  •  •  •    ^  ^(    0   ^^   'the 

order  statistics  of  a  sample  of  size  n  (>  k)  from  a  continuous 
distribution.   Let  us  denote 

H^^^^    ...  X^^^)   =   E{*(X-^  ...  Xj^)  I  0^3        (2.16) 

then  V(^)  <_  kV(4))/n  provided  that  V((t))  is  finite. 

Proof:   Since 

the  lemma  is  directly  derived  from  the  following! 

If  X-j ,  Xp,  ...   are  independently  and  identically 
distributed,  then 

cov  [H\    ...  X.,  X.^^  ...  X^),  .(X^  ...  X.,X^^^..,X^^_.)] 

li  v(<:>)  . 

In  order  to  prove  the  above  statement  we  can  assume 
with  no  loss  of  generality  that  ¥(*)  =  1. 
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Define 

c.    ^ov  {.(X,..  Xj  Xj^,  ..  X^)    *(Xi  ..  Xj  Xj^^,  ..  X^i^.j)} 

J  =   .  .  •  «  ^  K.  • 

h   =  k-j,  are  independently  and  identically  distributed 
random  variables. 

V  il^^i\   •••  ^y   Hi   •••  ^hi) 
h 

-  YZ  ^i^   •••  X-  X.  ,^  ...  x.,^)}  >  0 


i=l 
implies  that 


2h  +  2h(h-l)c.  -  2h^c  .  ,   >^  0 

Cj+i  1  kTJ  II  +  (k-j-l)cj)  ,  j  =0,1,2,... k-1. 


(2.17) 


Since  Cq  =  0,  it  is  proved  by  induction  that  c  .  j<  jAj 
J  =  l,...,k.  Q.E.D. 


From  the  lemma  we  have 


v(S?)  1  I  v(^,) 


hence 


nV(0?)  1  ^  +  e  .  (2.18) 

Thus  it  was  shown  that  0^  is  e-efficient  for  n  >  k  for 
all  f  in  ^  . 

The  coefficients  c   in  0     given  by  (2.^)  are 
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dependent  on  the  assiimed  density  t,    but  if  n  is  larger  than 
k  it  can  tie  also  estimated  from  the  sample.    Define 

^ap  =  ^(Y(a)  I  °n)        for  a  =  p 

=  ^°^  ^^(a)^^O)  I  °n)     "^^^ 

where  V(  ]  )  and  Gov  (  |  )  denote  the  conditional  variance 
and  covariance  given  the  order  statistics.   They  are  given  by 

s    =  E(Y?  ^  I  0  )  -  T^ 
aa     ^  (a)  '   n^     a 

P  Xf. s  -  T 
a   (i)     a 

(2.19) 


^ap  =  "^^(a)^YO)  I  °n)  '  ^a^p 


■^^  X/,.  \  X/^.  ^  -  T_  Tq   for  a  <  P 


^Pa  ~  ^ap 


waere 


pij    p  fY     X     Y    -X   J    i-lVl  ,1-i-l^P-a-l  n-j^k-p 

for  a  ;^  i,    p^n-j+k 
0  otherwise         (2.20) 

And  if  the  distribution  has  the  finite  fourth  order  moment, 
and  from  the  lemma  above. 
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Hence  for  fixed  k. 


^3  -^    °a3        i"  prcv-^-,..,,  (2.22) 


as  n  — >    T.      Hence   if  vi      ^ --'-•-- 

then  c  — >  c   in  probability,  and  if  we  put 


^*  =  5:  c   T 

■^ —  n        I 


a  a 


§*   and  0^  are  asymptotically  equivalent,  or  more  precisely. 


.S(r-§°)  =Z  (=„-<=;)  .^T^->0 


a 

in  probability.   Thus  ^*  is  asymptotically  equivalent  to 
an  e-efficient  estimator  9     for  any  f  in  ^. 

If  J-    is  properly  defined,  we  can  define  a  bounded 
open  subset  C  of  k;  dimensional  Euclidean  space  such  tnat 
for  any  f  in  3^  the  corresponding  vector  {c  )  lies  in  C. 
Then  we  can  define 


c  =  c     if   {c  )  €  C  ; 
a    a         a      ' 

otherwise  we  determine  c   so  as  to  minimize 

a 


H  51  s^  c   c 


ap  ^a  ^p 


under  the  condition  5"  c  :=  1  and  fc  )  e  C  where  C  denotes 
the  convex  closure  of  C. 
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Then  it  can  be  shown  that 

E(  ^  (c  -  c*)2)  _>  0  and  Y"  (c  -c*)^  is  bounded; 
therefore  if  we  define 

f*  =  Z  S„  T^  .  (2.25) 

nE{r-tO)2<E[{i:(?„-=a)')  '  ^  "  T^)  ] 

->  0 

as  n  — >  00  .  Hence 

nV(e^)  1  ^  +  e'  (2.2^) 

^f 

when  n  is  large  enough  for  any  e'  >  e. 

Thus  if  the  class  ^  satisfies  a  set  of  regularity 
conditions  thus  far  mentioned^  for  any  £  >  0  we  can  find 
a  k  for  which  (2.9)  is  satisfied  with  e/2,  for  all  f  in  ^, 
and  then  if  we  take  n  large  enough,  (2.24)  is  satisfied 
with  e'  =  £,  for  all  f  in  3^  . 

Finally,  if  we  choose  a  sequence  of  positive  numbers 
e  X  0,  and  construct  a  sequence  of  estimators  Q        in  the 
way  described  above  with  e  =  e   ,  then  we  have 

lim  n^{\)   ->    1  ,  (2.25) 

m— >oo  f 

where  {n  )  is  a  sequence  of  integers  tending  to  infinity. 


m 
Also  as  was  shown  by  Chernoff  et  al.  [8],  /FT  (S^.-  9 

is  asymptotically  normal   N(0  l/I^.)  when  k  — >  oo  ,  hence  if 

the  distributiors  in  ^  satisfy   the  conditions  mentioned 
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in    '^^     ir.ifonnly,    then  if  we   take   k   large   enough,   we   can 

obtain 

*  2 

sup   1p   [  y^TTT  (e°-e)    <  t)    -    f    ^-e-^   /2  dx|      <     e, 

for  arbitrary    e  >  0  for  all   f  in   \7^  . 

Then  again   if  we    take  n  large   enough,   we  have 

^  2 

sup    I  P^{    y^rr:   (^  -  0)    <   t)    -    f     -^  e"^    /^  dx|      <     2e   , 

-00    ' 

for  all  f  in  N?"  .   Hence  when  n  is  large  enough,  the  differ- 
ence between  the  distribution  of  9  and  that  of  the  normal 
becomes  smaller  than  any  prescribed  number,  uniformly  for 
all  the  distributions  in  ^  . 

Assuming  that  the  distribution  is  unknown,  the 
variance  of  d     is  also  unknown,  but  is  estimated  asymptoti- 
cally in  the  following  way. 

When  k  is  large. 


and 


But 


kl 


V(^?)  .  ^ 


V(fp  =  E(V(fj  0^))  +  V(E(0^|  0^)) 


Hence 


K^  E(V(fj  0^))  ^  V(^°)  (2.26) 
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,^     ,  *      * 


Since  V(0^  |  0^)  =  X  ^  ^ct  ^^   ^a^    '    ^^  ^^   replace  c^  by  c 


<^0- 


we  have  an  estimator  for  Y{9„)    given  by 


And  if  we  can  assume  asymptotic  normality,  we  have  the 
asymptotic  confidence  interval  for  9. 

We  shall  introduce  a  slight  modification  for  9 
Since  the  distribution  is  assumed  to  be  symmetric. 


°ap    \-a+l,k-p+l 


Hence  we  can  use  s  ^  =  (s  „  +  s        o.n)/2  instead  of  s  _, . 

ap    ^  ap    k-a+l,k-p+l'' '  ap 

^* 
We  shall  denote  by  9      the  estimator  of  the  form 

9*  =  y  c    T 

- —  a  a 

(2.28) 

where   c   =  ^I  ^^  /  Jl  ^  s^^  ,    Ib^'P)  =Is  „3-\ 
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§5.        'arlo  Results. 

We  performed  severu^  .  .i,.  ,-^  ;;v_.,.i._  jbi-^.^  uxperiiu 
for  •     distributions. 

1.   :.':^rmal 

50/6   Contaminated  normal 

f  (x)   =  0.95  X  -^  e"^   /^  +0.05     X  -^^-  e"^   ^^ 
)/5f  vTBtt 

lO/o  Contaminated  normal 

2  2 

f  (x)  =  0.90  X  -^-  e"^  /^  +  0.10  X  -^^-   e"^  ^^ 

4.  Logistic 

X 

F(x)  =  J  f(t)  dt 

-00 

5.  Double   exponential 


00  1  +  ^^P    (-  Tf  ^) 


f(x)    =  i-e-l^ 


6.  Tukey's    [l4]    special  distribution 

X  =  0^  [u-°-i  -  (i-u)-°-i] 

where  U  is   a   uniform  random  variable   over    [0  1]. 

7.  t-distribution  with  d.f.    2 

f(x)    = 


(l+x2)5/2 
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8.   Cauchy 


9.   Rectangular 


1  +  2x 


f(x)  =    -^        ,  lx|   <  y^ 

yT2 


10. 

Triangular 

f(x)  = 

/2 
/  3  " 

2 
"  3 

X  , 

X 

^1 

/I 

11. 

Quadratic 

1      I  I  2 
=  15-  (5-1^1  )  ^   1  i  1^1  1  3 

X  =  2(U  +  V+W)-3.   U;,V  and  W  are  independent 
uniform  random  variables. 

At  the  first  stage  we  computed  the  estimates  from 
random  samples  of  size  n  =  10,  from  each  of  the  populations 
listed  above  except  t(d.f.  2)  ,  for  k  =  3,-^,5,6,7,8.  N  =  1000 
samples  are  drawn.   Both  symmetrized  and  unsymmetrized 
estimates  are  computed  and  we  estimated  the  variances  of 
estimators  for  each  case.   The  results  are  listed  in  Table  1 
and  Table  2,  where  relative  efficiencies  are  calculated  by 

Ef f (9)    =      V3J"iaJ^ce  of  BLUE  —^ 
Variance  of  9 


1.   The  variances  and  coefficients  of  BLUE's  were  given  to  me 
by  Dr.  Morris  Meisner  who  obtained  the  numerical  results 
m  preparing  for  a  joint  paper  by  Birnbaum,  Laska  and 
Meisner  [6]. 
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The  results  s?    '  ;  priority  of  the 

symmetrized  est  .etrized  ones,  so  that 

in  ti.  experiments,  only  .^  -:■_ 

used. 

We  also  computed  the  means  of 


M  = 


a   3    ^ 
and  computed  the  ratio 

Mean  of  M 


Variance  of  9 


which  is  tabulated  in  Table  3.  The  values  are  fairly  uniform 
between  different  shapes  of  distributions,  except  for  Gauchy, 
and  are  uniformly  smaller  than  the  asymptotic  value  (n-k.)/k. 

As  a  rough  approximation  based  on  an  intuitive 
conjecture,  we  put  them  approximately  equal  to 


c  =  jn-kHn-k+l)  ,   . 


'k  ~    k(n-l) 
We  computed  studentized  values  using  this  approximation  by 

A 

e 


V^/Cj^ 


and  computed  cumulative  distributions  (Table  ^).   Again  the 
distributions  for  different  shapes  are  fairly  uniform,  also 
except  for  the  case  of  Cauchy,  in  which  the  distribution  is 
more  concentrated  than  the  other  cases.   And  the  distribu- 
tions (except  for  Cauchy)  of  the  studentized  values  are  seen 
to  be  fairly  close  to  t-distributions  with  degrees  of  freedom 

n-k+1. 
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Also  for  fixed  k  =  6,  we  computed  the  variances  and 
relative  efficiencies  of  estimators  for  n  =  10,15,20,  with 
N  =  1000,  the  results  being  given  in  Table  5. 

In  view  of  the  results  obtained  by  the  first  set  of 
Monte  Carlo  experiments,  it  seems  to  be  certain  that  the 
relative  efficiencies  of  our  estimators  are  high,  and  this 
implies  that  the  correlation  between  BLUE'S  and  our  estimator 
must  be  close  to  1.   This  implies  that  we  can  estimate  the 
difference 

v(^*)  -  v(Slue) 

much  more  accurately  thatn  V(9  )  itself  from  the  mean  of 

(3.2) 


A 


BLUE 


Since  we  have  exact  values  of  V(9p^-|-,,p,)  for  several  shapes 
of  distributions,  we  can  estimate  V(0  )  by 

V(e*)  =  Mean  of  (^^^-^blue)  +  ^^^BLUe)  *       ^^'^^ 

In  the  second  stage  of  the  Monte  Carlo  experiments, 
we  used  this  "difference"  estimator,  and  also  computed  a 
95  percent  confidence  interval  based  on  the  sample  variances 
of  the  difference  (3.2)  using  normal  approximation  for  the 
means.   The  results  are  given  in  Tables  6,  7  and  8  for 
n  =  10,  15,  20,  and  for  the  cases  where  the  BLUE  could  be 
obtained.   One  exceptional  case  is  rectangular  where  instead 

°^  ^BLUE  =  ^^min  +  ^max^^'^'  ^  "  sample  mean  is  used  (except 
for  the  case  k  =  6). 
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._.   jrror  f rv  .  .  rvals 

based  c-      above  appr   *•  *'  ns  of  the  student!" 
and  numbers  of  the  c;  2000  trials  in  which  the 

true  value  falls  outsid       .tervals  with  nominal 
confidence  levels  99,  93,  95*  90  percent  are  tabulated 
in  Tables  9,  10,  and  11. 

Finally,  for  k  =  5  and  k  =  4,  since  the  coefficients 
of  BLUE'S  are  all  nonnegative  except  for  extremely  irregular 
distributions,  and  surely  so  for  those  shapes  included  in 
our  list,  we  tried  the  estimator  with  nonnegative  coeffici- 
ents, that  is  we  put  the  condition  that  c   >  0.   The  results 

ex  — 

for  this  estimator  are  tabulated  in  Table  12  for  n  =  5  and  10, 
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Table  1.   Variances  and  Relative  Efficiencies  of  the 

Symmetrized  Estimators,   n  =  10     N  =  1000 


k 


7 


8 


BLUE 


Normal 

.1056 

.1076 



.1109 

.1137 

.1164 

1 
.1187 

.1000 

94% 

92°o 

90°/) 

87°/) 

86^0 

84% 

100% 

CN    .05 

.1286 

.1282.' 

.1285 

.1309 

.1339 

.135^ 

.1196 

93 

94 

94 

92 

90 

89 

100% 

CN    .10 

.1555 

.1540 

.1507 

.1524 

.1544 

.1563 

.1358 

37 

88 

90 

89 

88 

87 

100% 

Legist 

.  0964 

.  0979 

.0991 

.1016 

.1034 

.1045 

.0934 

97 

95 

94 

92 

91 

89 

100% 

Dl.    Exp. 

.1630 

.1631 

.1608 

.1626 

.1596 

.1633 

.1399 

86 

86 

87 

86 

88 

86 

100% 

Spec . 

.1053 

.1060 

.1054 

.1078 

.1087 

.1099 

.0957 

91 

90 

91 

89 

88 

87 

100% 

Cauchy 

1.1910 

.8392 

.2290 

.2220 

.2024 

.2146 

.1631 

14 

19 

72 

74 

81 

76 

100% 

Rect . 

.  0947 

.0905 

.0935 

.0953 

.1050 

.1095 

.0455 

48 

50 

48 

47 

43 

41 

100% 

Triang. 

.1048 

.1073 

.1126 

.1163 

.1222 

.1259 

Quad. 

.1056 

.1064 

.1108 

.1136 

.1177 

.1194 
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Non-Sy 


: .  J     ri  ■ 


Lzed   Estimators 


n 


.e 
N 


k  = 

7 

LUE 

Normal 

• 

• 

.337^ 
30 

0:333 
28 

.559^ 
28 

J  (.-    r  ^7 

.1 

CN    .-^^ 

■ 

.3534 

.3728 

.375 

.  -x^-. 

53 

32 

32 

'J<^ 

100% 

CN    .10 

.    '^^ZiL 

• 

■  r^l  -7 

.4082 

.4088 

.4078 

.1558 

^3 

:)[ 

'     . 

53 

55 

55 

ioo°/6 

Legist 

.2109 

.2460 

.2. 

.2701 

.2765 

.2749 

.095^ 

hk 

38 

;>^ 

55 

5^ 

54 

\Q0% 

Dl.    Exp. 

.2764 

.3066 

.3222 

.3318 

•  5582 

.5540 

.1599 

51 

46 

43 

42 

41 

42 

100% 

Spec. 

.20^7 

.2351 

.2492 

.2506 

.2623 

.2609 

.0957 

47 

41 

38 

58 

56 

57 

100% 

Cauchy 

'-^+352 

.35^1 

.3281 

.5025 

.5020 

.5105 

.1651 

37 

46 

50 

54 

54 

55 

100% 

Rect. 

.4829 

.5984 

.6586 

.6877 

.6912 

.6871 

.0455 

9 

8 

■   7 

7 

7 

7 

100% 

Triang. 

.3086 

.3602 

.3876 

.4022 

.4022 

.4050 

Quad. 

.2744 

.3184 

.3^43 

.5551 

.5507 

.5615 
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Table  5.   The  Ratios  of  the  Conditional  Variances 
the 
10 


and  the  Variance  of  Estimators 
P   _  (n-k)(n-k+l) 


n 


k(n-ir 


k 


Normal 

GN  .05 

CN  .10 

Logist 

Dl.  Exp, 

Special 

Gauchy 

Rect, 

Triang. 

Quad. 


3 


2.017 
2.051 
2.066 

2.097 
2.207 
2.122 
2.136 
2.178 
2.007 
2.073 


4 


1.076 
1.224 
1.240 
1.247 

1.317 
1.267 

1.675 
1.274 
1.170 
1.219 


5 


2.074     1.167 


-I 


0.670 
0.699 
0.717 
0.706 

0.735 
0.719 
1.205 
0.683 
0.641 
0.681 
0.667 


0.334 
0.4o4 
0.416 
0.4o4 
0.428 
0.412 
0.674 
0.372 
0.362 
0.390 
0.390 


7 


0.183 
0.192 

0.199 
0.196 

0.210 
0.200 
0.289 
0.161 
0.169 
0.185 
0.190 


8 


0.0767 
0.0809 
0.0847 
0.0847 
0.0905 
0.0869 
0.1195 

0.0633 
0.0697 
0.0794 
0.0833 
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k  =  3 


r 

.■',.  ■ 

CN.IO 

Lo 

DE 

Spec 

Cauchy 

Rect 

Trl 

Quad 

-4.0 

1 

1 

1 

1 

1 

1 

0 

3 

1 

1 

-3.8 

4 

1 

1 

1 

1 

1 

1 

4 

2 

( 

-3.6 

5 

1 

1 

3 

1 

1 

1 

4 

4 

;-; 

-3.4 

7 

3 

1 

3 

1 

3 

1 

5 

5 

9 

-3.2 

8 

4 

1 

3 

2 

4 

1 

5 

9 

10 

-3.0 

11 

7 

4 

3 

2 

4 

2 

7 

12 

13 

-2.8 

13 

10 

0 
"^^ 

8 

4 

7 

3 

12 

17 

15 

-2.6 

15 

18 

14 

10 

6 

9 

3 

17 

25 

18 

-2.4 

20 

22 

18 

15 

14 

13 

3 

20 

30 

24 

-2.2 

29 

29 

23 

22 

24 

23 

5 

23 

5$ 

28 

-2.0 

40 

36 

34 

32 

35 

37 

7 

31 

46 

36 

-1.8 

52 

48 

48 

48 

50 

51 

18 

42 

65 

50 

-1.6 

70 

65 

67 

73 

68 

74 

36 

67 

79 

69 

-1.4 

102 

97 

94 

112 

91 

103 

59 

103 

104 

94 

-1.2 

138 

136 

134 

140 

113 

136 

101 

127 

129 

123 

-1.0 

170 

170 

169 

187 

155 

173 

146 

171 

173 

164 

-0.8 

227 

230 

231 

232 

209 

225 

190 

215 

227 

207 

-0.6 

289 

301 

304 

501 

280 

287 

271 

271 

287 

274 

-0.4 

352 

364 

371 

367 

353 

357 

344 

545 

339 

359 

-0.2 

435 

445 

449 

429 

422 

431 

425 

4l8 

414 

446 

0.0 

508 

509 

520 

494 

491 

504 

492 

484 

477 

513 

0.2 

590 

588 

591 

568 

561 

572 

575 

576 

559 

595 

0.4 

648 

658 

653 

642 

633 

632 

659 

667 

637 

671 

0.6 

710 

714 

717 

717 

710 

701 

733 

740 

713 

730 

0.8 

771 

773 

773 

776 

762 

771 

789 

793 

782 

784 

1.0 

822 

828 

828 

824 

826 

817 

862 

847 

852 

832 

1.2 

862 

863 

866 

866 

871 

861 

901 

875 

863 

859 

1.4 

907 

901 

901 

898 

904 

888 

935 

913 

899 

889 

1.6 

933 

931 

931 

925 

926 

926 

960 

936 

924 

917 

1.8 

951 

952 

951 

948 

942 

954 

976 

949 

935 

937 

2.0 

960 

964 

965 

961 

964 

970 

982 

962 

949 

951 

2.2 

968 

970 

973 

977 

979 

982 

985 

972 

968 

960 

2.4 

978 

980 

981 

984 

986 

988 

991 

979 

973 

970 

2.6 

984 

985 

988 

990 

990 

990 

996 

986 

980 

979 

2.8 

988 

990 

991 

993 

995 

994 

996 

990 

986 

984 

3.0 

991 

992 

993 

994 

998 

997 

997 

992 

990 

990 

3.2 

994 

993 

994 

996 

999 

997 

1000 

994 

993 

993 

3.4 

996 

995 

994 

997 

1000 

997 

1000 

997 

996 

994 

3.6 

996 

996 

997 

998 

1000 

999 

1000 

997 

997 

995 

3.8 

998 

998 

998 

999 

1000 

999 

1000 

998 

998 

997 

4.0 

998 

998 

998 

999 

1000 

999 

1000 

998 

998 

998 

Table  4a.   CDF  of  studentlzed  estimator, 


t  d.f.    8.      tQ^^^5  =   3.355 


^0.99  "  ^'^^^       ^0.975 


^0.95     =  ^'^^^       '0.90  =  1-597 


'0.75 


2.306 
0.706 
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k  =   4 


N 

CN.05 

ON.  10 

Lo 

DE 

Spec 

Cauchy 

Rect 

Tri 

Quad 

-4.0 

2 

0 

1 

1 

1 

2 

1 

3 

1 

6 

-3.8 

5 

2 

1 

1 

1 

2 

1 

3 

3 

6 

-3.6 

6 

3 

1 

2 

1 

2 

1 

5 

5 

7 

-3.4 

8 

4 

1 

3 

2 

3 

1 

5 

6 

8 

-3.2 

9 

6 

3 

3 

2 

3 

1 

5 

11 

10 

-3.0 

10 

7 

4 

3 

2 

4 

1 

6 

13 

12 

-2.8 

14 

11 

7 

5 

3 

5 

2 

11 

17 

14 

-2.6 

16 

16 

11 

9 

8 

7 

3 

15 

22 

16 

-2.4 

18 

21 

15 

15 

12 

18 

4 

20 

30 

23 

-2.2 

27 

27 

23 

19 

24 

24 

5 

24 

37 

26 

-2.0 

38 

34 

33 

36 

30 

32 

8 

30 

45 

34 

-1.8 

49 

46 

42 

54 

48 

51 

16 

52 

59 

51 

-1.6 

67 

61 

63 

74 

68 

68 

30 

68 

76 

68 

-1.4 

98 

91 

87 

104 

87 

101 

51 

86 

101 

95 

-1.2 

130 

128 

125 

135 

112 

135 

76 

117 

128 

121 

-1.0 

167 

168 

163 

177 

145 

172 

134 

156 

168 

157 

-0.8 

221 

228 

225 

229 

204 

218 

176 

203 

224 

202 

-0.6 

289 

297 

296 

294 

268 

282 

254 

257 

286 

280 

-0.4 

354 

363 

374 

365 

341 

351 

337 

335 

339 

353 

-0.2 

428 

438 

439 

431 

422 

420 

413 

417 

406 

440 

0.0 

508 

513 

522 

499 

491 

502 

489 

489 

475 

516 

0.2 

584 

580 

581 

578 

567 

575 

579 

584 

563 

590 

0.4 

653 

660 

654 

646 

636 

633 

667 

675 

636 

658 

0.6 

711 

719 

721 

717 

715 

708 

738 

752 

717 

737 

0.8 

770 

772 

774 

782 

779 

772 

808 

809 

779 

787 

1.0 

827 

834 

833 

829 

836 

827 

867 

854 

833 

836 

1.2 

869 

869 

874 

863 

872 

867 

905 

881 

873 

848 

1.4 

920 

905 

906 

896 

910 

900 

942 

919 

900 

892 

1.6 

935 

933 

932 

929 

925 

928 

965 

937 

924 

915 

1.8 

951 

956 

955 

949 

949 

948 

977 

946 

939 

938 

2.0 

959 

964 

966 

965 

966 

970 

982 

964 

951 

954 

2.2 

968 

969 

972 

974 

977 

982 

987 

971 

964 

960 

2.4 

979 

981 

984 

981 

983 

989 

992 

979 

972 

974 

2.6 

983 

986 

989 

991 

990 

991 

994 

986 

978 

9S 

2.8 

987 

990 

991 

993 

995 

996 

997 

990 

985 

986 

3.0 

993 

993 

993 

996 

995 

997 

999 

992 

991 

990 

3.2 

995 

994 

994 

997 

998 

997 

999 

997 

991 

993 

3.4 

997 

996 

996 

997 

999 

998 

999 

997 

995 

994 

3.6 

998 

997 

997 

998 

999 

999 

999 

997 

996 

995 

3.8 

998 

998 

998 

998 

1000 

999 

1000 

997 

997 

996 

4.0 

998 

998 

999 

998 

1000 

999 

1000 

997 

999 

998 

t   d.f.    7 


Table  4b. 

CDF  of 

studentized  estimator. 

^0.995  " 

3.499 

'0.99  -  ^-998   to. 975  =  ^-'^5 

^0.95  " 

1.895 

'0.90  '  l-*^^  t^,j^     =  0.711 
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k  =  «^ 


1 

CN.IO 

Lo 

DE 

Spec 

Cauchy  Rect 

Tri 

Quad 

-4.C 

'   5 

3 

3 

3 

2 

1 

3 

3 

6 

-3.? 

4 

3 

3 

3 

4 

3 

5 

5 

7 

-3.6 

5 

3 

3 

4 

3 

6 

7 

9 

-3.4 

^* 

■_ 

5 

4 

4 

4 

3 

8 

0 

10 

-3.2 

9 

7 

f 

5 

k 

6 

3 

10 

12 

10 

-3.0 

12 

10 

7 

6 

7 

3 

13 

15 

14 

-2.8 

16 

14 

^t 

10 

7 

11 

4 

14 

19 

18 

-2.6 

18 

18 

15 

14 

11 

13 

6 

17 

22 

20 

-2.4 

20 

23 

18 

19 

^^ 

21 

10 

22 

31 

23 

-2.2 

?, 

33 

32 

25 

26 

33 

12 

54 

43 

30 

-2.0 

42 

40 

42 

3- 

40 

15 

41 

52 

34 

-1.8 

51 

50 

49 

59 

51 

57 

23 

59 

64 

58 

-1.6 

70 

65 

70 

72 

,'-  r 

79 

39 

71 

82 

73 

-1.4 

100 

92 

-9 

104 

104 

59 

87 

104 

100 

-1.2 

131 

124 

121 

143 

111 

135 

95 

118 

137 

129 

-1.0 

177 

173 

167 

179 

146 

182 

140 

158 

179 

164 

-0.8 

226 

228 

222 

230 

211 

222 

190 

202 

229 

204 

-0.6 

228 

297 

295 

287 

269 

284 

252 

264 

284 

279 

-0.4 

355 

359 

367 

364 

337 

358 

322 

327 

349 

347 

-0.2 

421 

431 

433 

431 

4l8 

419 

407 

415 

403 

435 

0.0 

506 

510 

515 

500 

497 

497 

493 

496 

479 

514 

0.2 

579 

580 

586 

581 

568 

571 

585 

587 

566 

5PR 

0.4 

649 

650 

645 

645 

639 

638 

609 

676 

634 

655 

0.6 

711 

712 

715 

715 

707 

717 

737 

756 

712 

728 

0.8 

765 

767 

772 

774 

788 

772 

810 

815 

771 

785 

1.0 

822 

825 

820 

819 

827 

820 

855 

850 

823 

829 

1.2 

863 

860 

864 

861 

875 

859 

896 

880 

863 

859 

1.4 

901 

896 

895 

892 

905 

898 

940 

910 

903 

887 

1.6 

928 

925 

931 

920 

924 

930 

957 

924 

915 

911 

1.8 

939 

948 

953 

943 

945 

946 

971 

939 

930 
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88.0 

CI 

+.0030 

-94.7 

+. 0033 

-94.0 

+. 0034 

-93.9 

M 

.9042 

18.0 

.6478 

25.2 

.2368 

68.9 

Cauchy 

V 
CI 

00 

00 

.2296 
+. 0239 

62.6 
-76.6 

M 

.0957 

47.5 

.0915 

49.7 

.0948 

47.9 

Rect. 

V 

.0008 

46.9 

.  0021 

48.6 

.0049 

46.4 

CI 

+.0012 

-.48.] 

+.0020 

-50.8 

+.0031 

-49.6 
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M 

.1107 

/-  - 
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.1181 

84.7 

N 

V 

.00^9 

t   . 

> 

.0071 

82.8 

CI 

+.0022 
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+  .  UUiiij 

'VOO  .    / 

+. 0027 

-^86. 7 

M 

.1258 

95.1 

.1311 

91.2 

.J  ■ 

-9.9 

CN    .05 

V 

.0070 

93.1 

.0081 

89.3 

.  00«; 

CI 

+. 0026 

-^97.1 

+.0028 

-93. 2 

+.0029 

'V     ■         , 

M 

.1^38 

94.4 

.1491 

91.1 

.1517 

CN    .10 

V 

.0088 

92.6 

.0105 

89.2 

.0111 

CI 

+. 0029 

^96  A 

+.0032 

-93.1 

+.0033 

'v^92.1 

M 

.1014 

92.2 

.10^43 

89.6 

.1050 

89.0 

Legist 

V 

.0036 

90.5 

.0049 

87.7 

.0060 

87.0 

CI 

+.  0019 

-93-9 

+.0022 

-91.5 

+. 0024 

~91.l 

M 

.1699 

32.3 

.1682 

83.2 

.1706 

81.7 

Dl .    Exp . 

V 

.0298 

79.8 

.0260 

80.7 

.0307 

79.0 

CI 

+.0055 

^85.1 

+.0051 

-J&^.Q 

+.0056 

'x.84.4 

M 

.1060 

90.3 

.1080 

88.7 

.1086 

88.1 

Special 

V 

.0058 

88.3 

.0071 

86.7 

.0088 

85.8 

CI 

+. 0024 

~92.4 

+ .  0026 

~90.8 

+.0030 

-^90. 6 

M 

.2393 

68.2 

.2873 

56.8 

.3^99 

46.'^ 

Cauchy 

V 

.2913 

63.6 

00 

00 

CI 

+.0171 

-73.^ 

M 

.0961 

47.8 

.1061 

42.8 

.1087 

41.8 

Rect. 

V 

.0217 

^5.5 

.0090 

41.7 

.0092 

40.7 

CI 

+.oo^»7 

~50.2 

+.0030 

'J^4.1 

+.0030 

'v^3..0 
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Table  7A.   Estimated  Variances  and  Efficiencies  of  the 
Estimators  by  the  "Difference"  Estimation 
n  =  15    M=Mean   Var=Variance   C.I. =95% 
Efficiency  is  given  in  percentages. 


k  - 

3 

] 

5 

N  = 

2000 

2000 

2000 

Var.         Eff. 

Var.         Eff. 

Var.         Eff. 

M 

.0678 

98.4 

.0692 

96.4 

.0708 

94.2 

N 

V 

.0002 

97.5 

.0005 

95.0 

.0008 

92.5 

CI 

+. 0006 
.  0818 

-99.3 

+.0010 

'^.97.8 

+.0013 

-96.0 

M 

96.1 

.0814 

96.5 

.0819 

96.0 

ON. .05 

V 

.0012 

9^.3 

.0012 

94.7 

.0016 

93.9 

CI 

+  .  0016 

'N.98.0 

+. 0016 

-93. 5 

+.0018 

'x.98.1 

M 

.0977 

90.7 

.0956 

92.7 

.0948 

93.5 

CN    .10 

V 

.0035 

88.3 

.0027 

90.5 

.0032 

91.1 

CI 

+. 0026 

-93.2 

+. 0023 

--95.0 

+.0025 

'%.96.0 

M 

.0637 

97.0 

.0643 

96.1 

.0651 

94.9 

Logist 

V 

.0007 

95.2 

.0008 

94.4 

.0009 

93.1 

CI 

+.0012 

^98.9 

+.0012 

-97.9 

+.0013 

--96. 9 

M 

.1118 

78.7 

.1097 

80.2 

.1077 

81.7 

Dl.    Exp. 

V 

.0100 

75.7 

.0094 

77.1 

.0088 

78.6 

CI 

+.0045 

~82.0 

+. 0044 

-83.5 

+  .0042 

-85.0 

M 

.512'^ 

17.8 

.2876 

31.7 

.1484 

61.4 

Cauchy 

V 

00 

00 

.1069 

55.7 

CI 

+.0150 

— 

-^^8.3 

M 

.0599 

36.9 

.0545 

40.6 

.0516 

42.8 

Rect. 

V 

.0002 

36.5 

.0004 

39.7 

.0012 

41.5 

CI 

+. 0006 

-37.3 

+. 0009 

^\.b 

+.0016 

-Jxk.^ 

-31- 


Figure  7B.        n  =   I5 
Eff. 
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V 
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.  004  C 

CI 
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~96.6 

+. 0019 

-^.94. 5 

+.0020 

~y4. , 

M 

.0651 

94.9 

.0665 

92.9 

.0681 

90.7 

Logist 

V 
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93.5 

.0014 

91.3 

.0018 

89.0 

CI 

+  .0010 

'N.96.4 
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'^>94.6 

+  .0013 

'^    '  .      . 

M 

.1079 

81.6 

.1071 

82.2 

.1079 
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Dl.    Exp. 

V 

.0099 

79.2 

.0105 

79.7 

.0117 

7'  c .  • 

CI 
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^4.0 

+.0032 

--84.7 

+.0034 

'V.84. 

M 

.1299 

70.1 

.1128 

80.8 
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. 

Cauohy 

V 

.0405 

66.9 

.0304 

77.1 

.0392 

78.2 

CI 

+.0064 

-75.8 

+.0055 

'^.85.0 

+. 0063 

-87.9 

M 

.0496 

44.6 

.0506 

43.5 

.0510 

43.3 

Rect. 

V 

.0058 
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.0022 

42.4 

.0028 

41.9 

CI 
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^e.Q 

+.0015 

^h.Q 
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Table  8A  .   Ectiiaated  Variances  and  Efficiencies  of 

the  Estimators  by  the  "Difference"  Estimation. 

n  =   20        M=Mean     Var-Variance      CI, 95%      Eff.    in   % 


k  = 

5 

4 

5 

N  = 

2000 

2000 

2000 

M 

Var.         Eff. 

Var.         Eff. 

Var.        Eff. 

.0506 

98.8 

.0513 

97.4 

.0521 

96.0 

N 

V 

.0001 

98.0 

.0002 

96.2 

.0003 

94.5 

CI 

+. 0004 

-99.6 

+.0006 

-^98. 6 

+.0008 

-97.5 

M 

.0610 

95.9 

.0606 

96.5 

.0606 

96.5 

ON    .05 

V 

.0008 

94.1 

.  0007 

94.6 

.0007 

94.6 



CI 

+  .0012 

-97.8 

+.0012 

'x.98.6 

+.0012 

-98.6 

M 

.0715 

92.0 

.0696 

94.5 

.0686 

95.3 

CN    .10 

V 

.0019 

39.5 

.0013 

92.3 

.0012 

93.7 

CI 

+.0020 

-^.94 . 6 

+. 0016 

-^96.7 

+. 0016 

'x.98.0 

M 

.0480 

96.3 

.0481 

96.1 

.0484 

95.5 

Logist 

V 

.0004 

94.6 

.0004 

94.4 

.0004 

93.7 

CI 

+.0009 

^93.1 

+.0009 

-97.9 

+. 0009 

-97.3 

M 

.0839 

75.9 

.0809 

78.7 

.0785 

31.1 

Dl .    Exp . 

V 

.0067 

71.6 

.0058 

75.6 

.0058 

77.8 

CI 

+. 0036 

-79.2 

+. 0034 

'x.82.2 

+. 0034 

'x.84.8 

M 
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29.8 

.1435 

43.7 

.0949 

66.2 

Cauchy 

V 

.0142 

62.6 

CI 

+. 0054 

-^70. 2 

M 

.0442 

29.2 

.0395 

32.9 

.0364 

35.7 

Rect . 

V 

.0001 

23.9 

.0003 

32.2 

.0006 

3-4-7 

CI 

+. 0004 

-29.7 

+.0008 

-33.6 

i 

+.0011 

'^.36.8 
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n  =   P'^ 
Eff.    ^: 


k  = 

6 

7 

8 

N  = 

4000 

4000 

4000 

Var.         Eff. 

Var.         Eff. 

Var.        Eff. 

M 

.0527 

94.9 

.0536 

93.4 

.0543 

92.1 

V 

.0005 

93.7 

.0006 

92.0 

.0007 

90.6 

CI 

+ . 0007 

~96.l 

+. 0008 

'x.94.8 

+. 0009 

-93.6 

M 

.0607 

96.4 

.0619 

94.5 

.0625 

93.6 

CN    .05 

V 

.0008 

95.0 

.0011 

93.0 

.0012 

92.0 

'  r 

+. 0009 

-97.8 

+  .0010 

-^96.1 

+. 0011 

-95.3 

M 

.0690 

95.4 

.0699 

94.1 

.0703 

93.6 

CN    .10 

V 

.0012 

95.9 

.0013 

92.7 

.0019 

91.8 

CI 

+.0011 

-97.0 

+.0011 

-95.6 

+.  0014 

-95.5 

M 

.0^83 

95.7 

.0499 

92.6 

.0504 

91-7 

Logist 

V 

.0005 

94.3 

.0007 

91.1 

.0008 

90.1 

CI 

+. 0007 

-97.1 

+. 0008 

'v.94.1 

+. 0009 

-93.3 

M 

.0783 

81.3 

.0768 

83.0 

.0768 

83.0 

Dl.    Exp. 

V 

.0056 

78.9 

.0055 

80.4 

.0057 

80.4 

CI 

+.0024 

-85.9 

+. 0024 

'v.85.6 

+. 0024 

--85. 6 

M 

.0857 

73.3 

.0802 

78.3 
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80.6 

Gauchy 

V 

.0115 

70.5 

.  0094 

75.4 

.0087 

77.7 

CI 

+.0054 

-76.5 

+. 0031 

-^l.H 

+.0029 

-83.7 

M 

.0556 

58.7 

.0554 

38.9 

.0326 

39.9 

Rect. 

V 

.0051 

36.7 

.0011 

37.8 

.0013 

38.6 

CI 

+.0018 

'x4o.9 

+.0010 

^0.1 

+. 0011 

-41.3 
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k  - 

level 

3 

4 

5 

6 

7 

8 

Normal 

99% 
98 

95 
90 

19 

40 

84 

189 

18 

30 

77 

171 

16 
32 

87 
187 

11 

89 

183 

20 

44 

112 

119 

61 
134 
237 

CN.05 

99  96 
98 
95 
90 

16 

31 

81 

188 

14 

28 

76 

156 

14 

28 

78 

168 

18 

27 

86 

172 

19 
40 

103 
196 

30 

56 

135 
225 

CN.IO 

99  96 

98 

95 

90 

12 

28 

75 
167 

11 

24 

57 

147 

13 

27 

74 

166 

18 

25 
81 

169 

17 

39 

102 

200 

30 
52 

134 

214 

Loglst 

99  ro 
98 

95 
90 

15 

28 

91 
190 

12 

20 

82 

167 

17 

30 

87 

177 

15 

31 

72 

168 

28 

52 

103 

199 

18 

36 
92 

202 

Dl.Exp. 

99  </o 
98 
95 
90 

5 

29 

79 

174 

6 
21 

66 
173 

11 
22 

83 
162 

9 

22 

63 
136 

17 

178 

15 

28 

82 

174 

Special 

99  96 
98 

95 
90 

15 

26 

84 

181 

12 

19 
72 

159 

17 

28 

76 

174 

15 

60 
167 

28 

49 

99 

190 

18 

35 

91 

198 

t  d.f.2 

99  7o 

98 

95 
90 

12 

20 

59 

143 

8 

15 

53 

125 

10 

63 
127 

12 

17 

58 

1^3 

25 
43 

88 
167 

15 

^? 

76 

174 

Cauchy 

99  96 
98 

95 

90 

3 

11 
18 
87 

4 

8 

32 

73 

6 
11 

47 

84 

8 
11 
33 
95 

17 

32 

65 

127 

12 
25 

136 

Rect 

99  9^ 
98 

95 
90 

28 

45 

98 

186 

23 
41 

80 
178 

31 
54 

109 
192 

31 
56 

96 
175 

44 

71 

131 

247 

29 

56 

132 

258 

Triang 

99  9^ 
98 
95 
90 

24 

52 

120 

228 

18 

44 

118 

197 

24 

39 
106 
208 

20 

28 

P 
163 

31 

60 
124 
226 

24 

50 

126 

243 

Quad 

99% 

98 

95 
-_2P 

27 

55 

121 

210 

26 
24 

107 
198 

29 
41 

109 
210 

19 

29 

89 

183 

35 

64 

136 
241 

23 

57 

125 

2?5 

Table  9.      Error  frequencies  of  confidence  intervals. 
n  =  10    N  =  2000 
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k  = 

level 

5 

4 

5 

u 

7 

8 

Normal 

00  OS 

95 
90 

20 
40 

103 
187 

18 

57 

92 

198 

18 

41 

96 

196 

15 

55 

^^7 

199 

22 

59 
94 

188 

^  r 

178 

00  96 
90 

15 
55 
83 

178 

14 
55 

85 
175 

16 

55 

85 

177 

14 

51 

87 

181 

19 

52 

50 

175 

18 

51 

66 

165 

CN-10 

99% 
98 

95 

90 

18 

55 

86 

175 

15 
55 

76 

165 

17 

57 

77 

160 

14 

50 

80 

175 

17 
55 

78 
169 

12 

24 

61 

156 

Legist 

99  96 

98 

95 

90 

18 

50 

78 

186 

16 

51 

74 

169 

15 

55 

75 

166 

14 

:>^ 

82 

154 

15 

67 

168 

17 

54 

75 

169 

Dl.Exp. 

99  9^ 
98 

95 
90 

14 

25 

71 

169 

10 

19 

67 

159 

9 
16 

69 
154 

10 

25 

58 
157 

ll 

65 

155 

10 

22 

62 

148 

Special 

99  96 
98 

95 
90 

14 
56 

180 

14 
52 

68 
164 

15 

52 

67 

162 

15 
55 

78 

146 

14 

29 

65 

160 

12 

29 

67 

154 

t  d.f.2 

9996 
98 

95 
90 

6 

29 

45 

151 

6 

20 

55 

155 

6 
21 

58 
151 

9 
26 

62 
118 

9 

20 

49 
155 

15 
22 

52 
124 

Cauchy 

9996 
98 
95 
90 

5 

9 

52 

85 

2 

6 

^° 
80 

2 

8 

55 

76 

6 
10 

54 
89 

2 
10 
50 
77 

4 

6 

29 

75 

Rect 

9996 
98 

95 

90 

20 

40 

78 

160 

17 

67 
157 

18 

67 

140 

20 

40 

82 

147 

19 

128 

26 

44 

79 

155 

Triang 

99  96 
98 

95 
90 

28 

46 

91 

191 

27 

59 

85 

182 

26 

40 

84 

180 

26 
41 

?^ 
165 

24 

58 

90 

177 

20 

54 

92 

176 

Quad 

99  9^ 
98 
95 
90 

28 

49 
106 
201 

28 

45 

105 

195 

25 

40 

104 

194 

20 
40 

99 

186 

24 

96 
186 

18 
54 

160 

Table  10.   Error  frequencies  of  confidence  intervals, 
n  =  15     N  =  2000 
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k  = 

level 

3 

4 

5 

6 

7 

8 

Normal 

99  yb 
98 

95 
90 

14 

37 

89 

192 

19 

35 

84 

186 

18 
32 

86 
188 

l4 

29 

95 
181 

20 

31 

85 

180 

19 

31 

89 

176 

CN.05 

990/) 
98 

95 

90 

13 

25 

85 

173 

11 
25 

85 
177 

13 
25 

89 
177 

15 
33 

84 

175 

17 
27 

83 
166 

17 
26 

87 
155 

CN.IO 

9996 
98 

95 
90 

13 
29 
80 

173 

11 

27 

78 

186 

10 

26 

84 
170 

15 

33 

85 

163 

80 
161 

13 

28 

80 

160 

Legist 

99  9b 
98 

95 
90 

15 

35 

92 

173 

12 

35 

167 

12 

29 

91 

168 

19 

^? 

76 

181 

12 

29 

82 
176 

12 

27 

73 

175 

Dl . Exp . 

999b 
98 

95 
90 

9 

17 

86 

184 

8 

18 

83 

175 

9 
19 
48 

170 

11 

29 

68 

147 

10 

19 

74 

165 

11 
20 

76 
154 

Special 

9996 

98 

95 

90 

15 
31 

85 
168 

11 

27 

89 
161 

11 

27 

88 
165 

17 

27 

76 

170 

10 
28 

77 
166 

10 
27 

82 
161 

t  d.f.2 

99  96 

98 

95 
90 

9 

17 

55 

139 

8 

20 

55 
129 

8 

21 

60 

129 

8 

18 

58 

134 

6 

22 

58 

119 

6 
21 

57 
131 

Cauchy 

999b 
98 

95 

90 

5 
10 

17 

57 

2 

7 
22 

58 

2 

8 
29 

67 

1 

5 

23 

77 

2 

5 

33 

67 

2 

5 
35 
73 

Rect 

99  9b 

98 

95 
90 

21 

3^ 

88 

154 

15 

31 

74 

144 

14 

30 

75 

135 

15 
24 

65 
135 

14 

27 

75 

121 

14 

31 

68 
124 

Triang 

99  9^ 
98 

95 
90 

28 

45 

97 

181 

25 
45 
87 

176 

24 

44 

90 

171 

19 

38 

84 

180 

23 

44 

90 

163 

21 

42 

92 
159 

Quad 

99  9b 
98 
95 
90 

23 

47 
109 
209 

22 

44 
108 
198 

23 

48 
114 
195 

16 

40 

93 

191 

20 

46 

100 

202 

20 

47 

107 

202 

Table  11.   Error  frequencies  of  confidence  intervals, 
n  =  20    N  =  2000 
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w 


fficiencies  of  Estimators 
ifficients. 


n  = 

5 

5 

10 

4 

.2121 
9^.3 

.2177 
91.9 

.2000 
100% 

.1016 
98.4 

.1055 
96. 

CN   .05 

.2482 
100 

.2527 
98.8 

.2497 
100 

.1223 
97.6 

97.6 

CN   .10 

.2975 
97.0 

.5024 
95.^ 

.2885 
100 

.1467 
92.6 

.1458 
95.2 

Logist 

.2005 
95.2 

.2088 
91.5 

.1910 
100 

.0965 
96.9 

.0978 
95.^ 

Dl.    Exp. 

.5^96 
90.6 

.5641 
87.0 

.5169 
100 

.1707 
81.9 

.1705 
82.0 

Spec. 

.212^4 
95.1 

.2218 
91.1 

.2021 
100 

.1041 
91.9 

.1050 
91.0 

Cauchy 

2.400 
50.9 

2.551 
48.2 

1.221 
100 

.9451 
17.2 

.7255 
22.2 

Rect. 

.2447 

58.5 

.2452 
58.2 

.1428 
100 

.0960 

47.8 

.0918 
49.5 
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§^ ,   Observations  on  Monte  Carlo  Results. 

The  following  observations  are  made  on  Monte  Carlo 
results  stated  above. 

1.  Relative  efficiencies  with  respect  to  BLUE  are  quite 
high  for  small  sample  sizes  n  =  10,  1^,    20,  and  even  for  n  =5. 

2.  Except  for  the  double  exponential,  Cauchy  and  rectan- 
gular, all  of  which  have  some  irregularities  to  be  discussed 
in  the  following,  the  relative  efficienties  are  uniformly 
higher  than  SO   percent,  and  quite  often  higher  than  95  /o. 
5.     Double  exponential  distribution  has  a  discontinuity 

in  the  derivative  of  the  density,  hence  one  of  the  regularity 
conditions  is  not  satisfied.   Although  the  BLUE  is  asymptoti- 
cally equivalent  to  the  median  and  is  fully  efficient  in  this 
case,  whether  0      is  asymptotically  fully  efficient  is  open 
to  question.   The  answer  might  be  affirmative,  but  it  cannot 
be  derived  from  the  combinations  of  general  statements 
mentioned  in  §2. 

4.     The  Cauchy  distribution  has  no  moments,  and  for 
sample  size  3  or  4,  no  unbiased  estimator  with  finite 
variance  exists.   This  seems  to  give  an  intuitively  natural 

explanation  for  the  low  efficiencies  for  the  cases  k  =  3 

.  /s^     ^  2 

and  4.   Also  in  such  cases  the  sample  variances  of  9   -9 

are  so  big  and  unstable,  that  they  may  well  suggest  that 
E{9*    )  is  actually  infinity?  hence  we  denote  the  facts  by 
putting  00  in  the  table.   Ev^n  for  the  case  k  _>  5,  it  is  a 
problem  whether  9      actually  converges  in  probability  to  9„, 
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because  the  nonexistence  of  the  moments  makes  the  general 
arguments  inapplicable. 

~.     For  rectangular  distribution,  the  variance  of  the 

-2  -1 

BLUE  (=  (X    +  X  .  )/2)   is  proportional  to  n   not  to  n 

as  is  usually  the  case,  hence  it  is  expected  that 

hence  the  relative  efficiency  will  be  approximately  k/h. 

Relative  efficiencies  with  various  k  values  depend 
on  the  shapes  of  the  distribution.   For  the  normal  and  also 
for  sufficiently  regular  distributions  such  as  the  5  percent 
contaminated  normal  and  the  logistic,  smaller  k  gives  higher 
relative  efficiencies.   For  longer-tailed  distributions 
larger  k  is  better,  although  too  large  k  gives  uniformly 
poor  efficiency. 

7.  The  pattern  of  relative  efficiencies  with  different 
values  of  k  for  each  distribution  is  fairly  stable  as  n  varies 
from  10  to  20.   For  most  cases  the  relative  efficiency  for 
fixed  k  increases,  but  not  always,  and  even  when  it  does, 

the  increase  is  not  steady  or  rapid. 

8.  The  gain  in  the  efficiency  by  symmetrization  is 
surprising  when  n  =  10.   It  may  decrease  as  n  increases,  but 
it  is  not  yet  investigated.   On  the  other  hand  the  effect  of 
restricting  to  non-negative  coefficients  is  negligible  when 
n  =  10,  which  fact  can  be  seen  by  comparing  Table  6  and 
Table  12. 
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9.  Confidence  Intervals  based  on  the  "studentization" 
described  above  show  fairly  uniform  error  probabilities 
except  for  Gauchy  and  the  t-distribution  with  d.f .  2. 
Also,  the  error  frequencies  are  nearly  uniformly  smaller 
than  the  prescribed  level,  which  implies  that  a  better 
approximation  to  the  distribution  should  still  be  found. 

10.  If  we  take  the  normal  case  as  the  standard,  confi- 
dence intervals  tend  to  be  biased  more  or  less  on  the 
conservative  side  when  the  distribution  is  long-tailed, 

and  seem  to  be  slightly  biased  in  the  opposite  direction  when 
it  is  short-tailed. 

Some  further  remarks  are  in  order. 

11.  The  variances  of  BLUE  may  not  be  a  good  criterion 
for  calculating  the  relative  efficiencies.   The  most 
proper  bound  will  be  the  variance  of  Pitman's  estimator 
which  gives  the  sharp  lower  bound  for  all  location  invariant 
estimators,  but  the  numerical  values  are  not  available. 
The  Monte  Carlo  study  now  under  way  by  V.  Mik^  (personal 
communication)  indicates  that  the  variance  of  the  BLUE'S 

is  not  very  far  from  that  of  Pitman,  and  generally  much 
closer  to  his  than  the  Cram6r-Rao  bound,  when  the 
discrepancies  are  large  (e.g.  for  such  cases  as  the 
double  exponential  and  the  Cauchy).   Hence  the  relative 
efficiences  measured  by  the  variance  ratios  with  BLUE 
may  be  a  little  higher  than  more  plausible  ones^  we  can 
expect  them  to  be  not  too  far  from  them  and  to  be  a 
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reasonable  approximation. 

12.  The  relative  effic      3  themselves  are  sometimes 
misleading,  because  they  should  be  Judged  also  in  reference 
to  those  of  other  simpler  estimators.   For  example,  relative 
efficiencies  of  the  sample  mean  (which  is  actually  e qual  to 
the  symmetrized  estimator  with  k  =  2)  should  be  taken  into 
account.  Table  13  gives  the  efficiencies  of  the  sample  mean. 
Reference  should  also  be  made  to  the  small  sample  efficiencies 
of  the  trimmed  mean  (Gastwirth  and  Cohen  [12]). 

13.  The  relative  efficiencies  could  also  be  compared  with 
the  maximum  efficiencies  of  linear  estimators  (Bimbaum 

et  al.  [6])  or  of  the  maximum  (invariant)  estimators 
(Mik^  [21]).   The  numerical  results  obtained  by  them 
suggest  that  there  remains  little  room  for  improving  the 
efficiencies  of  our  estimator  over  many  shapes  of  distri- 
butions.  For  example  when  n  =  20,  the  maximum  efficiencies 
of  linear  estimators  given  in  [6]  are 

Cauchy  versus  Normal  83.5  /^ 

Double  exponential  vs.  Normal  91.5  /^ 
10%  contaminated  normal  vs.  Normal  96.6  /6 
Cauchy  vs.  Rectangular  22.5  /o 

Normal  vs.  Rectangular  67.5  /o 

Those  values  are  not  much  better  than  the  values  attained 
by  our  method  of  a  favorable  choice  of  k. 
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Table  I5.   Relative  Efficiencies  of  the  Sample  Mean.  {%) 


n  = 

5 

10 

15 

20 

Normal 

100.0 

100.0 

100.0 

100.0 

CN   .05 

88.1 

85.4 

84.2 

83.6 

GN    .10 

80.1 

75.4 

73.9 

73.7 

Legist 

95.5 

93.5 

92.7 

92.3 

Dl .    Exp . 

79.2 

70.0 

66.0 

63.7 

Special 

85.7 

81.1 

* 

■x- 

Cauchy 

0 

0 

0 

0 

Rect. 

71.4 

45.5 

33.0 

2.60 
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lU.  choice  of  k  :  also  on 

the  possible  shapes  of  distributions,  but  roughly  the 
following  will  be  a  plausible  range, 

k  =  5 

10  5 

15  ^   •' 

20  5-7 

More  roughly,  k  =  5  seems  to  be  a  reasonable  choice  for 
the  range  of  sample  sizes  n  =  10  ~  20,  and  for  relatively 
regular  distributions  such  as  the  normal,  contaminated 
normal,  or  the  logistic,  we  can  guarantee  relative 
efficiencies  of  approximately 

93  %     for    n  -  10 

94  15 

95  20 
or  slightly  higher. 

15.    The  relative  efficiencies  of  the  confidence  intervals 
based  on  studentization  will  be  lower  than  those  of  estimators 
because  of  the  small  degrees  of  freedom  of  the  approximate 
t-distribution,  especially  when  n-k  Is  small.   But  the 
exact  relative  efficiencies  are  hard  to  assess,  because 
little  has  been  known  about  the  small  sample  shortest 
confidence  intervals  for  the  median  of  non-normal  distribu- 
tions.  Further  investigation  of  this  problem  is  yet  to 
be  done. 
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16.  The  regularity  conditions  required  to  establish 
the  asymptotic  full  efficiency  of  our  estimator  are  not 
explicitly  stated.   Although  it  is  quite  possible  to 
state  them  rigorously  by  putting  the  results  of  several 
previous  papers  together,  they  may  be  unnecessarily 
restrictive,  and  it  seems  to  be  desirable  to  find  a 
possible  unified  set  of  regularity  conditions  which  are 
not  too  restrictive. 

17.  The  same  idea  as  proposed  here  can  be  applied  to 
obtain  robust  procedures  in  other  problems  about  location 
parameters.   Application  to  the  testing  about  one-sample 
and  two-sample  location  parameters  is  discussed  in  another 
paper  by  the  author  [25].   Also,  the  estimated  coefficients 
of  BLUE  can  be  used  to  test  the  hypothesis  about  the  shapes 
of  distribution  (e.g.  normality)  against  symmetric  alterna- 
tives, which  will  be  discussed  in  another  paper  [26]. 

18.  After  this  work  was  completed,  the  author  obtained 
a  new  result  which  gives  a  theoretical  foundation,  and 
also  some  generalization,  to  the  method  proposed  here,  and 
which  will  be  given  in  a  paper  [27]  now  in  preparation. 
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